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Abstract 

Introducing the generalized, non-extensive statistics proposed by 
Tsallis into the standard s-wave pairing BCS theory of supercon- 
ductivity in 2D yields a reasonable description of many of the main 
properties of high temperature superconductors, provided some al- 
lowance is made for non-phonon mediated interactions. 

The discovery of superconductivity in the copper oxides in 1986 |^ and 
the subsequent race for even higher critical temperatures (125K by 1993 
[P, ^) raised hopes for the application of superconducting phenomenon at 
operating temperatures approaching room temperature. The inability of the 
Bardeen-Cooper-Schrieffer (BCS) model [^] to describe superconductivity 
in these materials satisfactorily, appears to indicate that we are dealing 
with a completely different class of superconductors. Various theoretical 
models have been proposed to explain this phenomenom, ranging from d- 
wave superconductivity ^ to models incorporating non-phononic coupling 
mechanisms [§, 10, 11 1. Different degrees of success have been achieved in 



explaining specific aspects of these high- Tc materials, but no inclusive model 
exists. 

A number of characteristics must be incorporated in any such model. 
Certainly the main common denominator in all high-Tc materials is the 
large anisotropy in the crystal structure, resulting in conduction electron 
states in the Cu02 planes being very nearly two-dimensional in character. 
It is therefore reasonable to assume that this is essential for the high critical 
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temperatures and that arguably the 'ideal' high- Tc superconductor is purely 
two-dimensional in character. It is, however, well known that in BCS, a 
pure 2D model presents problems as the phase loses its coherence due to 
fluctuations [^]. This calls for at least a quasi-2D model in all realistic 
cases. Although the results of measurements on the superconducting gap 
are somewhat varied, there seems to be convergence towards a ratio of -r^S- 
equal to somewhere between 6 and 8 which is not consistent with the uni- 
versal value of ~ 3.5 obtained for normal BCS superconductors. This is 



particularly so in the case of tunneling measurements [03, 14, 15, Oq], where 



ratios as high as 8.9 have been measured in Lai.85Sro.i5Cu04 [17|, but ap- 
pears to contradict optical measurements ||l^, |l^, 20 1 which seem to yield 
results closer to the BCS ratio for weak coupled superconductors. In spite of 
the apparent BCS behaviour of the tunneling currents, the violation of uni- 
versality is difficult to understand. Thirdly, experimental data suggests that 
the electronic contribution to the specific heat in cuprates does not exhibit 
the exponential form of normal weak coupled superconductors. It is rather 
of the linear form 7T, where 7 is of the order of a few mJ/mol-K^. Magnetic 
properties of the high- Tc superconductors also differ appreciably from their 
normal counterparts, most notably in their very high upper critical fields. 
These fields are extrapolated to be of the order of tens or hundreds of Teslas 
and are too large to be attained with other present technologies. Further- 
more it should also be noted that all high-Tc superconductors are type II 
superconductors. One of the strongest arguments for the phonon coupling 
scheme of the BCS model is the existence of an isotope eff'ect. In cuprate 
materials, where there is a strong doping dependence of the isotope effect, 
this requires some additional consideration. In optimally doped materials 
the effect is often strongly suppressed (Tc may scale as a power a = 0.1 of 
the mass), whereas in overdoped or underdoped materials it might be more 
prevalent, with the scaling in some cases even exceeding a = 0.5 |21]. 

Recently in statistical mechanics there has been interest in a generaliza- 
tion of the Boltzmann-Gibbs (BG) statistics to a non-extensive form pro- 
posed by Tsallis|jl| in which the former is recovered in an appropriate limit. 
This formalism has had considerable success in providing an appropriate 
mathematical framework for dealing with physical systems with long-range 
interactions. This is markedly so in the so-called stellar "polytropes" |22] 
where the usual Boltzman distribution functions yield unphysical results. A 
variety of other applications have also been considered along with the gen- 
eralization of many well-known theorems and principles, see for example the 
references in [23|. 

In this paper we propose and motivate the possibility that the cuprate 
oxide materials have an inherent underlying non-BG like character respon- 
sible for their high-critical temperatures and other unusual properties. We 
suggest that the s-wave BCS model is essentially correct in employing an 
effective weak coupling Hamiltonian that includes the kinetic energy of free 
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electrons along with a constant attractive potential between electrons of 
equal momentum and opposite spin j^. The BCS ansatz in a purely 2D 
form is used for the ground state, with a generalized form of the Fermi- 
distribution function |24] appearing at finite temperatures. We show else- 
where that such a 2D model is justified as the fluctuation-dissipation the- 
orem is no longer valid in the Tsallis formalism |25|. The assumption that 
the coupling mechanism is purely phonon mediated, needs to be modified 
to obtain reasonable fits to the experimental data, particularly in the case 
of very high temperature superconductng materials. 

The relevant question is, of course, why one would resort to changing 
the entropic measure in a high-Tc superconductor, and if one does, why the 
Tsallis formulation should be appropriate. A very compelling argument can 
be found from a consideration of the electronic specific heat. The more 
exotic coupling methods mentioned above, certainly may have the effect of 
removing problems indirectly related to the BCS-Hamiltonian such as the 
lattice instabilities predicted by Migdal |26]. The experimental evidence of 
an energy gap, however, requires an energy spectrum for the total energy of 
excitations, of the form 

El = el + A^ (1) 

where £ is the excitation energy and A the energy gap, independent of any 
proposed model. Assuming a Fermi distribution, the electronic specific heat 
capacity (C) at temperatures kT <^ A, can shown to have the following 
form 1^] 

_^ 

C - (2) 

One is therefore obliged to contend with an exponential form of the heat 
capacity in any model that relies on the Fermi distribution function used 
in BG statistical mechanics. One way to circumvent this is to introduce a 
different distribution function. 

It should be remembered that the most basic interaction in an electron 
gas is that of the Coulomb repulsion. This interaction is of infinite range and 
may also introduce correlation effects. One example of this is the condensa- 
tion of a true electron gas, of low enough density, into a Wigner lattice |25]. 
Very early attempts at explaining superconductivity, quite intuitively, but 
unsuccesfully used the Coulomb potential as a starting point [^, More 
recently models like those of Hubbard describing contributions due to 
repulsion between electrons sharing an atomic orbital, have been applied to 
superconductivty. The condensation of electrons into its superfluid state is 
ultimately the result of an effective, attractive interaction. The extent of 
the effect of Coulomb repulsion might not be obvious and in spite of the fact 
that we do not explicitly take it into account in the effective interaction, we 
argue that it is precisely the effects of this long-range interaction which may 
render the system more suitable for description by generalized statistics. 
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The generalized entropy postulated by Tsallis ^ takes the form: 

Sg = k^^MPi^lH g 3fJ) (3) 

where w is the total number of microstates in the system and p is the 
associated probabilities. It has been shown that this entropy obeys the 
usual properties of concavity, equiprobability, positivity and irreversibility 
and Shannon additivity and preserves the Legendre transformation stucture 



of thermodynamics |32|. It is also straight forward to verify that the usual 



BG entropy is recovered in the limit q ^ 1. Associated with equation (H) is 



the generalized Fermi distribution given by |24] 



/. = 1 • (4) 

[l + /3(g-l)efc]— + 1 

Once again the usual Fermi distribution of BG statistics is recovered in the 
limit q ^ 1. It is this generalized distribution, we will henceforth assume, 
that excited independent quasi-particles obey at finite temperature. In 2D 
the BCS gap equation for weak coupled superconductors at finite tempera- 
ture 1^ is given as 

\ = -Yl 3kp-^i'^ - fkl - fki) (5) 

k ^ 

where g is the coupling constant, Ej. is the quasi-particle energy given by 
equation (IT]) and fka is the Fermi distribution function. Replacing f^a by 
equation (^) and the sum by an integral over the density of states we obtain 
in 2D a generalized form for the gap equation: 

1 _ /- de [l+/?(g-l)V£2 + A2]^ -1 
N{d)g J Vi^TA^[i+/3(g_i)Vi2^A^]i^ +1 

where N(0) represents the density of states. For q=l in 3D, one recovers 
the gap quation in the standard BCS theory and, of course, a description 
of normal superconductors in the weak coupling limit. As required, equa- 
tion @ is independent of the distribution function and at zero temperature 
reduces to 

-J—=f ^£ (7) 

where Aq is the zero temperature gap. 

Let us, as a starting point, choose as cutoff to the gap equation the 
Debye frequency Hcod- The analytical solution to equation (0) is 



Ao = {hiUD + Jhiol + A2)e W . (8) 
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Using experimental values for the gap and Debye frequency, one can solve 
for N{0)g at zero temperature. The problem then reduces to finding that 
q in the gap equation, which yields a vanishing gap at the critical temper- 
ature. With q thus fixed, the temperature dependence of the gap can be 
determined. 

Consider again the gap in equation (P) at Tc, independent of the par- 
ticular choice of statistics. A change in the integration variable e — > eksTc 
may always be made such that the integration limits are from to 
Integrating by parts yields: 

^ = (l-2/(|))l„|-/*(.-2F(.))l„. * (9) 

where F(e) is the indefinite integral of f(e). Taking the integral on the right 
to the left hand side and dividing by (1 — 2/(^)) yields some number ip 
which depends, of course, on the choice of /(e) . Exponentiating both sides 
yields: 



n = One'' (10) 

which clearly demonstrates that the isotope effect is preserved in the BCS 
formulation irrespective of the particular form of statistical mechanics em- 
ployed. 

It is at this stage appropriate to note that despite reasonable physical 
arguments in the case of normal superconductors, the choice of the Debye 
frequency as a cutoff is mathematically quite arbitrary. The temperature 
dependent solution to the gap equation in fact converges to a well defined 
value for any cutoff of sufficient magnitude. It is only when the cutoff is of 
a magnitude of the order of the gap itself, that deviations occur. In normal 
superconductors the Debye frequency is greater than the gap by ~ 10^ and 
thus adequate. This is, however, certainly not the case in high-Tc super- 
conductors where the gap at zero temperature may even be larger than the 



Debye frequency, e.g. in Tl2Ba2Ca2Cu30io ||3^ (see for Debye frequen- 
cies of other high Tc superconductors). Abandoning the Debye frequency 
as a cutoff is tantamount to acknowledging that the electron-phonon inter- 
action is not the only interaction involved in forming the condensed state. 
Changing the cutoff will, of course, influence the isotope effect. It might be 
argued that the observed suppression in the isotope effect may be a conse- 
quence of this. 

It is interesting to note that the convergence of the gap at greater cutoffs, 
is accompanied by a simultaneous convergence of q. Consider a transforma- 
tion of the gap in equation (^) via the substitution e = ^bTcc ^ -^^^ 
define the cutoff in terms of multiples of the energy gap, e.g. nAg. Then 
for a given value of the gap to critical temperature ratio, = m, one can 
show that the gap equation reduces to a form dependent only on the ratio 
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-^^j^, the number n and q, and not on either the gap or critical temperature 
directly. 

1 ^ r^^ de' [1 + {q - l)e' /2]^^ - 1 ^^^^ 
N (0)5 io e' [1 + (g _ iy/2] + 1 

Thus, if for a specific ratio of Aq to Tc, the temperature dependent gap is 
convergent independent of cutoff, then q must have a uniquely defined value. 
It therefore seems appropriate to define the cutoff in terms of the conver- 
gence of q as this can be specified independent of the critical temperature 
pertaining to a specific material. 




Figure 1: The ratio -jTpfr vs. q is given by the solid curve. For comparison 
the dashed curve representing ^^^^ vs. q is included. 

A relevant question now is whether a generalization of the BCS univer- 
sality condition of ^ 3.52 exists for the generalized statistics of Tsallis. 

Clearly the generalization must be q dependent because of the dependence 
of Tc on q and reduce to the BCS universality condition for q=l. In Fig.l 
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a graph of versus q is given. This ratio does not deviate apprecia- 

bly from 3.5 which suggests the following generalization of the universality 
condition 

4^ - 3.52 (12) 

A more detailed analysis might lead to replacing q by some function of q 
which should be ~ q. 




Figure 2: The normalized electronic specific heats y^— vs. ^ for various 

choices of -jz^j^, q and cutoffs where in each case Cmax corresponds to the 
maximum value of C. The corresponding values of the normalized gaps are 
shown in the inset. 

The electronic specific heat capacity may be expressed as 

C = T% (13) 

where the entropy is given by 
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Figure 3: Comparison of the temperature dependence of experimental 
values of A from Ref lH for Bi2Sr2CaCu208 with = 85K and ^ = 6.2 
with the theoretical results for g=1.71 



S, = 2ks[ ^''^f^f -/(^^''^ + - ^(^j) - - f^^'^^'h (14) 

and we use A^(0) = for the density of states in 2D and the bare electron 
mass TJie in all cases. The effect of the cutoff on the specific heat is shown in 
Figure 2 for LSCO (Tc = 36K). The linear nature of the specific heat can 
be seen over most of the superconducting region with huo = 390K as the 
cutoff. Clearly the situation deteriotes as the critical temperature is reached. 
Increasing the cutoff to TOAq yields a more linear form of the specific heat 
near Tc with little change in the shape of the gap (see inset in Fig. 2). 
Changing the cutoff from htor) to 70Ao increases the 7 from 5.5 mJ/mol-K'^ 
to 10.7 mJ/mol-K^ with a slight change in q (from 1.752 to 1.679). Note 
that in spite of the fact that the slopes of the specific heat may be altered 
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by using an effective mass (7 ^7); tliese values are not in disagreement 
with many of the experimental results which seem to lie between 3 and 12 
mJ/mol-K^ for most cuprates [p4| For comparison the ratios of ^^j, = 4.0 
and 8.0, the corresponding results for the specific heats and the gaps are 
shown. Note the specific heat is much more linear when larger cutoffs are 
used and that the nonlinearity around Tc disappears almost completely for 
1^ = 8.0. In this case 7 = 20.0 J/mol-K^ . 

In Fig. 3 the experimental data obtained by Briceno and Zettl |13] for 



the temperature dependence of the gap of Bi2Sr2CaCu208 with Tc = 85K 

9 A 

and -^^j^ = 6.2 are compared with the theoretical results with g=1.717. In 
spite of a linear specific heat (7 = 13.9 mJ/mol-K^ using mg) which is quite 
consistent with experiment jsj], the shape of the gap differs slightly from 
that obtained experimentally. 

Due to the type II nature of the cuprate oxides, a full analysis of the crit- 
ical magnetic field would require including contributions to the free energy 
resulting from the flux lattice as well as the mutual energy of the surface 
current and flux line currents. We therefore present only a qualitative anal- 
ysis of the thermodynamic critical field, and make no attempt at solving 
either the lower- (i?ci) or the upper- (-Bc2) critical fields. 

In 2D, the thermodynamic critical field is given by 

^ = Fn- Fs (15) 

where and Fg are the free energies in the normal- and superconducting 
states respectively and //q is the permeability of free space. The free energy 
in the superconducting state is given by 



Fs = N{m^?[[l + - 1] - 2iV(0) / de{^^-^)f,{m (16) 

AkT r 

- ^^(0) J ^^[(1 - ^'?(^^))' - (1 - ^'?(^^))]- (1^) 

Fn is obtained by setting A = and using the value of q determined from 
solution of the corresponding gap equation. This implies that in our gen- 
eralized BCS model the normal state is also described by the generalized 
statistics of Tsallis since Be must vanish at Tq- Figure 4 shows the ther- 
modynamic field obtained for a superconductor with §^jr = 6.0 and Tc = 
36K . The inset shows the behaviour of the free energies. They differ from 
the weak coupled normal superconductor case where the free energy is al- 
ways negative except at T = where it vanishes. The difference, however, 
remains a positive definite quantitiy and no problems arise in calculating 
Be- The zero temperature result of figure 4 agrees exceptionally well with 



an extrapolation based on experimental results due to |35]. They find, for 
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a Lai.846Sro.i54Cu04 crystal with Tc = 35K that Be = 0.251 T, while our 
result is 0.234 T. 

The Josephson quasi-particle current was first treated in detail by Shapiro 



et al. |^6[ (the other terms are considered elsewhere |37|). The quasi-particle 
current between two different superconductors is given by the integral: 

1 /■°°_, , „, .. 0(|cj| - Ai)0(|(j + ey| - Aa) 



bRn J-oo (u;2 _ Af) 2 [(u; + eVy - A^] 2 

(18) 

where V is the applied voltage, e the electron charge, Rat the junction resis- 
tance and Al, A2 the energy gap on either side of the tunneling junction. In 
Figure 5 we show the result of evaluating equation (|l8|) at 5K for a junction 
with identical superconductors having = 6.0 and using the generalized 
distribution function with g=1.679. Of interest is the fact that Iqp is not neg- 
ligible below an applied voltage of 2A/e. The equivalent calculation using 
BG statistics yields a vanishing contribution in this region. Included in the 
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Figure 5: ^vs.|^ (inset) obtained from the quasi-particle contribution to 
the Josephson tunnehng current in the region eV< 2A, as shown in the main 
figure. 

figure is the derivative ^ . Although we do not obtain all of the structure ob- 
served experimentally in LSC0||3^, the existence of a Igp contribution raises 
the possibility that interference due to other admixtures, perhaps d-wave or 
Andreev bound states, will produce the sought after structure. Analysis 
of the other contributions using generalized statistics does not appreciably 
alter the results in the eV < 2A region. The use of generalized statistics, 
in the BCS equations, however, allows the experimentally observed discon- 
tinuity in the quasi-particle tunneling current to be interpreted as being 
coincident with the energy gap at low T, whilst consistently predicting the 
gap to vanish at the correct critical temperature. 

In conclusion we have shown that a simple 2D s-wave BCS pairing model 
which incorporates generalized statistics can provide an adequate description 
of many of the main features of high temperature superconductivity. Non- 
phonon mediated interactions, however, appear to play a role as the critical 
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temperature increases. 



References 

[1] C. Tsallis, Journal of Statistical Physics 52(1/2), 479 (1988). 

[2] J.G. Bednorz and K.A. Muller, Z. Phys. B 64, 198 (1986). 

[3] S.N. Putilin, E.V. Antipov, O. Chmaissem, and M. Marezio, Nature 
362, 226 (1993). 

[4] A. Schilling, M. Cantoni, J.D. Guo, and H.R. Ott, Nature 363, 56 
(1993). 

[5] J. Bardeen, L.N. Cooper, and J.R. SchriefFer, Physical Review 108(5), 
1175 (1 Dec. 1957). 

[6] J.F. Annett, Cont. Phys. 36(6), 423 (1995). 

[7] P. Muzikar, J. Phys.: Condens. Matter 9, 1159 (1997). 

[8] A.S. Alexandrov and J. Ranniger, Physical Review B 23, 1796 (1981). 

[9] CM. Varma, Physical Review Letters 63, 1996 (1989). 

[10] P.W. Anderson, Science 235, 1196 (1987). 

[11] W.A. Little, Physical Review A 134, 1416 (1964). 

[12] P.C. Hohenbcrg, Physical Review 158(2), 383 (10 Jun. 1967). 

[13] G. Briceno and A. Zettl, Sol. State. Comm. 70(11), 1055 (1989). 

[14] H.L. Edwards, J.T. Markert, and A.L. de Lozanne, Phys. Rev. Lett. 
69(20), 2967 (16 Nov. 1992). 

[15] Z. Zhang and CM. Lieber, Phys. Rev. B 47(6), 3423 (1 Feb. 1993). 

[16] G.T. Jeong, J.L Kye, S.H. Chun, S. Lee, S.L Lee, and Z.G. Khim, Phys. 
Rev.. B 49(21), 15 416 (1 Jun. 1994). 

[17] T. Ekino, T. Doukan, H. Fujii, F. Nakamura, S. Sakita, M. Kodama, 
and T. Fujita, Physica C 263(1-4), 249 (1996). 

[18] D.A. Bonn, J.E. Greedan, CV. Stager, T. Timusk, and M.G. Doss, 
Phys. Rev. B 35(16), 8843 (1 Jun. 1987). 

[19] L. Degiorgi, E. Kaldis, and P. Wachter, Sol. State Comms. 64(6), 873 
(1987). 



12 



[20] R.T. Collins, Z. Schlesinger, R.H. Koch, R.B. Laibowitz, and T.S Plas- 
kett, Phys. Rev. Lett. 59(6), 704 (10 Aug. 1987). 

[21] J.R. Waldram, Superconductivity of Metals and Cuprates (lOP Pub- 
lishing Ltd, 1996). 

[22] A.R. Plastino and A. Plastino, Physics Letters A 174, 384 (1993). 

[23] A.R. Plastino, A. Plastino, and C. Tsallis, J. Phys. A: Math. Gen. 27, 
5707 (1994). 

[24] F. Buyukkilic, D. Demirhan, and A. Gulec, Physics Letters A 197, 209 
(1995). 

[25] H. Uys and H. G. Miller, Preprint, University of Pretoria (2001). 
[26] A.B. Migdal, Sov. Phys. - JETP 7, 996 (1958). 

[27] E.M. Lifshitz and L.P. Pitaeskii, Statistical Physics, Part 2, vol. 9 (Perg- 
amon Press, 1980). 

[28] E.P. Wigner, Phys. Rev. 46, 1002 (1934). 

[29] W. Heisenberg, Z. Naturforsch. 2a, 185 (1947). 

[30] H. Koppe, Z. Physik 148, 135 (1957). 

[31] J. Hubbard, Proc. Roy. Soc. (London) A276, 238 (1963). 

[32] E.M.F. Curado and C. Tsallis, J. Phys. A. 24, L69 (1991). 

[33] A. Hoffmann, P. Lemmens, G. Gutherodt, V. Thomas, and K. Winzer, 
Physica C 235-240, 1897 (1994). 

[34] D.M. Ginsberg, Physical properties of high temperature superconductors 
II (World Scientific, 1990). 

[35] Q. Li, M. Suenaga, T. Kimura, and K. Kishio, Phys. Rev. B 47(17), 11 
384 (1 May 1993). 

[36] S. Shapiro, P.H. Smith, J. Nicol, J.L. Miles, and P.F. Strong, IBM J. 
Res. Dev. 6, 34 (Jan. 1962). 

[37] R.E. Harris, Phys. Rev. B 10(1), 84 (1 Jun. 1974). 

[38] L. Alff, A. Beck, R. Gross, A. Marx, S Kleefisch, T. Bauch, H. Sato, 
M. Naito, and G. Koren, Phys. Rev. B 58(17), 11197 (1 Nov. 1998). 



13 



List of Figures 

1 The ratio t^t^ vs- q is given by the solid curve. For compar- 



ison the dashed curve representing -;^r^T\ vs- q is included 



The normalized electronic specific heats -^r^ vs. ^ for var- 



ious choices of q and cutoffs where in each case C, 



corresponds to the maximum value of C. The corresponding 



values of the normalized gaps are shown in the inset. 



Comparison of the temperature dependence of experimental 



2Ao 



values of A from Ref [13 1 for Bi2Sr2CaCu208 with Tc = 85K 





4 The critical magnetic field Be vs. T obtained from free ener- 


gies given in the inset. 




5 ^vs.f^ (inset) obtained from the quasi-particle contribution 


to the Josephson tunneling current in the region eV< 2A, as 
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